Field dependence of the Spin State and Spectroscopic Modes of Multiferroic BiFeOs 
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The spectroscopic modes of BiFeOs provide detailed information about the very small anisotropy 
and Dzyaloshinskii-Moriya (DM) interactions responsible for the long- wavelength, distorted cycloid 
that appears below Tn = 640 K. A microscopic model that includes two DM interactions and easy- 
axis anisotropy is able to closely describe both the zero-field spectroscopic modes as well as their 
splitting and evolution in a magnetic field. While only six modes are optically active in zero field, 
all of the zone-center modes are activated by a magnetic field. The close agreement with experiment 
suggests that the proposed model provides the foundation for future technological applications of 
this multiferroic material. 

PACS numbers: 75.25.-j, 75.30.Ds, 78.30.-j, 75.50.Ee 



Due to the coupHng between electric and magnetic 
properties, mutliferroic materials are among the most im- 
portant yet discovered. With a multiferroic material used 
as a storage medium, information can be written electri- 
cally and then read magnetically without Joule heatingi. 
Hence, applications of a room-temperature multiferroic 
would radically transform the magnetic storage industry. 
Because it is the only known room-temperature multifer- 
roic, BiFeOs continues to attract intense interest. 

Although its ferroelectric transition temperatur^ 
Tc ~ 1100 K is much higher than its Neel transition 
temperature^— Tn ^ 640 K, the appearance of a long- 
wavelength cycloid^i^ — with a period of 62 nm enhances 
the ferroelectric polarization below Tn- The induced po- 
larization has been used to switch between magnetic do- 
mains with an applied electric field^i^i^. 

Progress in understanding the microscopic interactions 
in BiFeOs has been greatly accelerated by the recent 
availability of single crystals for both elastic and inelastic 
neutron-scattering measurements. By fitting the spin- 
wave (SW) frequencies above a few meV, recent inelas- 
tic measurements^^' ^^ have determined the antiferromag- 
netic (AF) nearest-neighbor and next-nearest neighbor 
exchanges Ji ^ —4.5 meV and J 2 = —0.2 meV, indi- 
cated in the pseudo-cubic unit cell of Fig. 1(a) with lattice 
constant a ~ 3.96 A. In the presence of strair^ii^ or a 
magnetic fieldi^ii^ above He ^ 20 T, those exchange in- 
teractions produce a G-type AF with ferromagnetic (FM) 
alignment of the S = 5/2 Fe^+ spins within each hexag- 
onal plane. 

The cycloidal magnetic order in BiFeOs below He 
is produced by the much smaller anisotropy and 
Dzyaloshinskii-Moriya (DM) interactions. Because the 
wavelength a/\^S of BiFeOs is so large, inelastic 
neutron-scattering measurements do not have the q- 
space resolution needed to resolve the cycloidal satellites 
at q = (27r/a)(0.5 ± (5, 0.5, 0.5 T ^)^ on either side of the 
AF wavevector Qo = (7r/a)(l, 1,1). Below 5 meV, inelas- 
tic measurements at Qo reveal four broad peaks, each of 
which can be roughly assigned to one or more of the SW 
branches averaged over the first Brillouin zoneiiii^. 

By contrast, THz spectroscop y^^! ^^ provides very pre- 
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FIG. 1: (a) The pseudo-cubic cell with exchange interactions 
Ji and J2 as well as the polarization direction z' cutting 
through two hexagonal planes, (b) For any of the three mag- 
netic domains, a schematic of the spins in zero field show- 
ing their rotation about y^ Due to the DM interaction 
D^ = D'z' , spins rotate by r about z^ in the x'y' plane, (c) 
The wavevector parameter S versus field for domain 1 (solid) 
and domains 2 and 3 (dashed) with vertical dashed lines show- 
ing their critical fields, (d) The energies of domains 1 (solid) 
and domains 2 and 3 (dashed) versus field. Both (c) and (d) 
assume that the field is along [0,0,1]. 



cise measurements for the optically-active SW modes 
at the cycloidal wavevector Q. For polarization along 
z' = [1,1,1] (all unit vectors are assumed normal- 
ized to 1), the three magnetic domains have ordering 
wavevectors Qi = (27r/a)(0.5 + (5,0.5 — (5,0.5) (domain 
1), Q2 = (27r/a)(0.5 + (5,0.5,0.5 - (5) (domain 2), and 
Qs = (27r/a)(0.5,0.5 + (5,0.5 - 6) (domain 3). The local 
coordinate system {x^^y^^z'} for each domain is indicated 
in Fig.l. 

In zero field, the four observed mode frequencies were 
recently described by a modeP^ with two DM interac- 



tions and easy-axis anisotropy K along z^ While the 
DM interaction D along y' is responsible for the cycloidal 
period, the DM interaction^"— D' along z^ produces a 
small tilt^ r in the cyloidal plane, shown in Fig. 1(b), 
which alternates in sign from one hexagonal plane to the 
next. In the canted AF phase above Hc^ D' produces a 
weak FM momenUiii^ of about Mq = 0.03/iBy'. 

This Letter shows that the microscopic model de- 
scribed above also predicts the mode splitting and evolu- 
tion of the spectroscopic modes with field. Due to mode 
mixing, all of the SWs are optically active in a magnetic 
field. The mode splitting and field dependence allows 
us to unambiguously assign the spectroscopic modes of 
BiFeOa. The close agreement between predictions and 
measurements suggests that the proposed model can pro- 
vide the foundation for future work on BiFeOs. 

In a magnetic field H = Hm along m, the spin state 
and SW excitations of BiFeOs are evaluated from the 
Hamiltonian 

H = -Ji ^ Si • S, - J2 ^ S, • S, - i^^(S, • z'f 



ihj) 



ihj)' 



-L> V / • (Si X S.) 

■'^^Rj=Ri+ax, 



^'Er r, (-1)^"'/^z'.(S,xS,) 



(1) 



While the nearest- and next-nearest neighbor exchange 
interactions Ji = —4.5 meV and J2 = —0.2 meV can 
be obtained from the SW dispersion between 5.5 meV 
and 72 meV^^-^^, the small interactions D, D\ and K 
that control the cycloid can be obtained from THz spec- 
troscopy below 5.5 meV (44.3 cm~^). 

For a given set of interaction parameters and magnetic 
field, the spin state of BiFeOs is obtained by minimizing 
the energy E = (H) over a set of variational parameters. 
With the same spin states in hexagonal layers n and n+2, 
we parameterize the spin states in layers n = 1 and 2 as 



^^/(R) = A(^)(R)sin/icosr(^) sin(27r(5i?^//a + 7r0 

+ sopi':\ (2 

Sy^iK) = A(^)(R)sin/isinr(^^ sin(27r(5i?^//a + 7^''^) 

+ SOPy^\ 

5,.(R) = A(")(R)cos/iF(")(R) + 8opi?\ 



(2) 



(3) 



where 



(4) 

F^^)(R) = ^C2i-i cos(2(2/ - l)7r6R^^/a) 

+ Yl ^21 co8{4l7rSR,^/a + F^")) (5) 

1=1 



and we take Ci = 1. Notice that the unit vectors p^^^ 
and the tilt angles r^^^ can be different for layers 1 and 



2. Four phases 7}^^ and 72 parameterize S'a;/(R) and 
Sy'CR). In zero field, the higher odd harmonics C2i-\-i>i 
in F^^\Rx') are produced by either the anisotropy K 
or the DM interaction D' . In addition, even harmonics 
C21 are produced by the magnetic field. Because the 
coefficients Ci fall off rapidly with /, we neglect harmonics 
above 1 = 4. The additional variable F^^^ on each layer 
allows the even and odd harmonics to be out of phase. 
On layer n and site R, the amplitude A'^^^IV) is fixed 
by the condition that |S(R)| = S^ which is satisfied by a 
quadratic equation for A^^\lV). The lower root is used 
for layer 1; the upper root is used for layer 2. 

Fixing 5 = l/g, where g ^ 1 is an integer, E is 
minimized over the 17 variational parameters (/i, r^^\ 
Ji \ p(n)^ pM ^ C/<4, and sq) on a unit cell with q sites 
along x' and two hexagonal layers. Another minimization 
loop is then performed over q to determine the cycloidal 
wave vector as a function of field. In zero field, q = 222. 
We verify that the corresponding spin state provides at 
least a metastable minimum of the energy by checking 
that the classical forces on each spin vanish. Another 
check is that the SW frequencies calculated below are all 
real. 

We emphasize that the variational parameters are not 
free but rather are functions of the interaction parameters 
D^ D' ^ and i^, and the magnetic field H. In zero field, 
the spin state reduces to the one used in Ref.[19]. A 
much simpler variational form for the spin state would 
have been possible were the field oriented along the high- 
symmetry axis z' = [1,1,1]. The complex form given 
by Eqs.(j2]|4|) is required because the experimental field 
is oriented along [0,0,1]. Above Hc^ the canted phase 
with (5 = can be more simply parameterized by taking 
Q>i = 0, r(^) = r(2), r(^) = 0, and p(^) = m, thereby 
reducing the number of variational parameters to 7. 

Taking m = [0, 0, 1], |m • x'| and |m • y'| are the same 
for domains 2 and 3. Therefore, their equilibrium and 
dynamical properties are identical. In Fig. 1(c), we plot 
5 versus field for the three domains. Notice that the 
critical field He = 16.4 T for domains 2 and 3 is lower 
than i^c = 20.2 T for domain 1. The variation of He with 
m was recently reported in Ref. [15| and was predicted by 
Le Bras et al?"^ for a purely harmonic cycloid. In zero 
field, all three domains have the same energy. But in a 
nonzero field, domain 1 has lower energy than domains 2 
and 3, as seen in Fig. 1(d). At 5 T, the energy difference 
between domains is about 0.9 /ieV/site. 

Generally, the spin-spin correlation function ^'^/^(q, a;) 
may be expanded in a series of delta functions at each 
SW frequency C(j^(q): 



S'c/3(q,cj) = ^^(^ 



.(q))^i^^(q), 



(6) 



which assumes that the SWs are not damped. The mode 
frequencies ujm{^ and the corresponding intensities Sv^^ 
are solved by using the 1/ S formalism outlined in Ref.|25| 
and in Appendix A of Ref.[19| for non-collinear spins. 
With (5 = 1/g, the unit cell contains 2g sublattices. 
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FIG. 2: The mode frequencies versus q for (a) T and (b) 6.9 
T in domain 1. Optically- active modes at 77 = ^ are denoted 
by filled circles, inactive ones by white circles. 



Only some of the SW modes are optically active 
with non-zero magnetic resonance (MR) matrix elements 
I ((5 1 Mq, 1 0)1, where M^ = 2/iB "^Zi ^ioi is the magnetization 
operator, |0) is the ground state with no SWs, and \d) 
is an excited state with a single SW mode at the cy- 
cloidal wavevector. A subset of the MR modes are also 
electromagnetically (EM) active with nonzero matrix el- 
ements |((5|P^"^|0)|, where P^"^ is the operator for the 
electric polarization induced by the cycloid. Expressions 
for |(^|MJO)| and \{^\P'^'^\^)\ are provided in Appendix 
BofRef.[l|. 

For zero field with 8 = 1/222, we adjusted^^ the inter- 
action parameters of BiFeOs to fit the observed spectro- 
scopic modesii at 17.5, 20.1, 21.5, and 27.3 cm~^. Using 
the weak FM moment Mq = 0.03/iB above He to fix2S 
D' = MqJi/ ji^S^ we obtained the interaction parame- 
ters D = 0.107 meV, D' = 0.054 meV, and K = 0.0035 
meV. We now employ those interaction parameters to 
describe the field dependence of the spectroscopic modes 
in BiFeOs. 

In Fig. 2, we plot the SW frequencies versus q = 
(27r/a)(0.5 + 7^,0.5 - 7^,0.5) for domain 1 and ilf = 
or 6.9 T. To label the spectroscopic modes at q = Q or 
7^ = (^, we have modified the notation of de Sousa and 
Moore^i, who studied the case where D' and K vanish 
and the cycloid is purely harmonic. In an extended zone 
scheme, they labeled the SW modes at wavevector nQ 
as ^n and ^n- Corresponding to excitations within the 
cycloidal plane, ^n = *^i|^| is a linear function oin. The 
out-of-plane modes satisfy the relation ^^ = <^i\/l + n^- 
Due to the higher harmonics of the cycloi d^^i^^ produced 
by D' or K^ ^n>o and ^n>o each split into two modes 
that we label as ^i^n and ^^i^q- 



Only <l>o, ^0, ^1 



n>0 

(1) 



^f\ 
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^i\ and ^2 ^^^ optically 



active in zero field. Those six modes are indicated by the 
filled circles in Fig. 2 (a). Although the lowest-frequency 
$0 niode has an enormous SW intensity and MR matrix 
element |(^|M2y/|0)|, it is not detectable by THz mea- 
surements. At 21.5 cm-\ ^^^^ with \{5\M^>\{))\ > 
is also EM active. The nearly-degenerate ^0 and ^i ^ 
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FIG. 3: The measured spectroscopic frequencies with rf field 
hi (circles) or h2 (triangles). The predicted mode frequencies 
from domain 1 (solid) and domains 2 and 3 (dashed) are also 
shown. We conjecture that contributions from domains 2 and 
3 stop at 6 T, indicated by a dashed vertical line. 



modes with significant |(^|Ma^/|0)| and |((5|M2y/|0)| matrix 
elements, respectively, are responsible for the observed 
spectroscopic peak^^ at 17.5 cm~^. While ^^ i^ acti- 
vated by the 3Q harmonic of the cycloid, which mixes 
^f^ with '^[^\ ^0 and <l>^^^ are activated by the tilt of 
the cycloid out of the x^z^ plane, which mixes ^0 with 
^^^^ and ^^^^ with $0. 

In nonzero field, all of the zone-center SW modes are 
optically active with nonzero MR matrix elements. Many 
of these modes also have nonzero EM matrix elements. 
The MR- active modes below 5.5 meV are indicated in 
Fig.2(b). Notice that ^0 and ^\ are split by the mag- 
netic field. 

With m = [0, 0, 1], the measured field-dependenc o^^i^^ 
of the spectroscopic modes is plotted in Fig. 3. The 
frequency-dependent magnetic field was aligned either 
along hi = [1,-1,0] (circles) or h2 = [1,1,0] (trian- 
gles). These r/ fields couple to the MR matrix elements 
|((5|M • h^|0)|. The observed transition to the canted 
phase occurs just below 20 T. Due to instrumental lim- 
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FIG. 4: The evolution of the matrix elements |(^|Ma|0)|//iB 



with field for (a) ^\ and (b) ^j^ Mn domain 1, where ol = x 
(solid), y' (medium dash), or z' (small dash). On the right 
are the intensities Sx'x\^^ and Syfyf{6) (thin solid curves) of 
those modes. 



itations, no THz data is available for fields above 12 T 
and frequencies below about 12 cm~^. 

This figure also contains the predicted spectroscopic 
frequencies for domain 1 (solid), and domains 2 and 3 
(dashed). We believe that the energy difference between 
domains is responsible for depopulating domains 2 and 3 
above about 6 T, indicated by a dashed vertical line. To 
reflect this behavior, we have cut off the predicted mode 
frequencies of domains 2 and 3 above 6 T. 

The agreement between the measured and predicted 
mode frequencies in Fig. 3 is remarkable. In agreement 
with predictions, ^o and ^[ ^ are slightly lower in domain 
1 than in domains 2 and 3. The predicted splitting of 
^2 is clearly seen in Fig. 3. Also in agreement with 
predictions, ^^^ Ms slightly lower in domains 2 and 3 than 
in domain 1. The only feature that remains unexplained 
by our model is the field-independent mode observed at 



about 16.5 cm~^ midway between ^i ^ and ^o- Notice 
that spectroscopic modes never cross with field due to 
their coupling. Because the constant 16.5 cm~^ mode 
appears to be immune to mode repulsion, we speculate 
that it may be an optical phonon. 

The curves in Figs. 3 include several modes that be- 
come optically active only in a magnetic field H. At 10.8 
cm~^, ^[ ^ has no active MR matrix elements in zero 
field but develops significant matrix elements |(J|Ma^/|0)| 
and \{6\My^\0)\ in a field, as shown in Fig.4(a). The SW 

intensity Sx'x'i^) of ^i grows like H^. Experimentally, 
^l appears above about 3 T. 

The selection rules governing the observation of the 
spectroscopic modes also change with external magnetic 
field H. At 20.4 cm-\ ^f^ has \{S\M^^\0)\ > at zero 

field. Since hi -z' = but h2 -z' ^ 0, ^^ ^ is only observ- 
able in r/ field h2, given by the triangles in Fig. 3. But 
Fig. 4(b) shows that |((^|Maj/|0)| increases linearly with 

("2) 

field, which explains why ^-^ appears in r/ field hi when 

H > 0. The intensity Sy'y'{5) of ^\ ^ decreases only 
slightly with ilf up to 6 T. 

Recall that we use the same interaction parameters 
previously obtained for zero field^^. Therefore, it is not 
surprising that the quantitative agreement with measure- 
ments, although quite good, is not perfect. In particu- 
lar, the measured spectroscopic frequency of the canted 
phase above He is about 4 cm~^ lower than predicted. 
Improvements might be acheived by using slightly higher 
values^^ of L)^ 

To conclude, the close agreement between the pre- 
dicted and observed spectroscopic frequencies with field 
leaves no doubt about the origin of those modes. This 
agreement suggests that the present model, with DM 
interactions along y' and z' and easy-axis anisotropy 
along z', can provide the foundation for future studies 
on BiFeOs and may lay the groundwork for its eventual 
technological applications. 
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